A quantum beam splitter for atoms. 
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An interferometric method is proposed to controUably split an atomic condensate in two spatial 
components with strongly reduced population fluctuations. All steps in our proposal are in current 
use in cold atom laboratories, and we show with a theoretical calculation that our proposal is very 
robust against imperfections of the interferometer. 

PACS numbers: 3.75.Fi, 39.20.+q 
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I. INTRODUCTION 

Cold atoms are routinely applied in spectroscopy and 
in atomic clocks, and matter wave interferometers have 
been used to probe with unprecedented precision the ef- 
fects caused by acceleration, gravity or interaction with 
atoms or electromagnetic fields ||]. The experimental 
resolution in these experiments is linked to the precision 
with which the relevant signal can be measured and to 
the intrinsic uncertainty of the quantity measured. For 
example, in Ramsey interferometry applied in atomic 
clocks, one tunes a radiation frequency to obtain a 50 
% population of both states of the clock transition, but 
if the atoms are uncorrelated, in addition to the measure- 
ment errors one will inevitably observe differences scaling 
as ^/N between the populations of the two states, where 
N is the total number of atoms (binomial distribution) . 
So-called spin squeezed states have been shown to ex- 
ist, for which it is possible to significantly reduce this 
uncertainty without going to exceedingly large numbers 
of atoms ||. A number of proposals has been made 
to produce such spin squeezed states by interaction with 
squeezed light ^, |^], by quantum non-demolition mea- 
surement of the atomic populations pi m , and by light- 
induced ^ or coUisionally induced |P,^o|r]ll[| interaction. 

In this paper we propose a method to realize a beam 
splitter for a one-component condensate, causing a split- 
ting of the condensate in two spatially separated compo- 
nents with a better matching of occupancies than in the 
binomial distribution, resulting from a normal splitting 
of the particles. Such splitting has indeed been proposed 
to occur if one adiabatically raises a potential barrier in- 
side a single condensate, and experiments have confirmed 
the effects of coUisional interactions on the population 
statistics of a condensate split by a periodic potential 
[ p^ . This dynamics is governed by the time scale for 
which the system is able to adiabatically follow the low- 
est energy state due to the collisional repulsion among 
atoms, and this time scale may be very lon g, m aking an 
experimental implementation very difficult |13| . 

In contrast, we propose a fast method working in four 
steps: (i) apply a normal splitting to create a binomial 
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FIG. 1: Atom-interferometer realized via Bragg scattering of 
the condensate. 



distribution of atoms separated in space, (ii) make use 
of the collisional interaction in each spatial component 
to cause in few msecs a non-linear phase evolution of the 
different amplitudes, (in) reflect the atoms so they again 
overlap in space, and (iv) remix them on a second beam 
splitter, so that the resulting separate components are 
populated by a sub-binomial distribution. 

In Sec. II, we present a theoretical analysis of the in- 
terferometer applied to a Bose-Einstein condensate. In 
Sec. Ill, we discuss the results ideally obtained with our 
proposal, and in Sec. IV, we analyze the effect of a mis- 
match in overlap in the interferometer. 



II. BRAGG INTERFEROMETER 

The principle of the proposal is sketched in Fig. ^, 
where the vertical arrows indicate Bragg diffracting laser 
fields with appropriate detunings and phases. We sug- 
gest to use Bragg diffraction to split and to recombine 
the atomic clouds because this method has already been 
succesfuUy demonstrated in experiments. 

Our proposal combines ideas for light squeezing in a 
Kerr-nonlinear interferometer fl^ with ideas for spin 
squeezing of two-component condensates using the col- 
lisional interactions among the atoms In a sin- 
gle condensate, the collisional interaction leads to an ef- 
fective quadratic term in the Hamiltonian proportional 
to N{N — 1), where N is the total number of atoms. 
That term affects only the global phase of the con- 
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densate, but if the condensate is divided in two com- 
ponents with no mutual interaction, it is replaced by 
''T-i{ni — 1) + n2{n2 — 1), having a significant effect on 
a superposition of states with different populations ni 
and n2 in the two components. The lack of interaction 
between the components could be due to a tuning of the 
collisional properties of the atoms, or, as suggested in 

due to a spatial separation of the atoms. In 
we carried out simulations with the so-called positive-P 
method which retains the full many-body character of 
the problem. Our calculations were well accounted for 
by a simple two-mode description, and we shall therefore 
apply such a description in the following analysis. 



A. Internal and external dynamics 

The atoms are initially in a single component BEG, 
i.e., all atoms populate the same one-particle wave func- 
tion ilj{r). Counterpropagating laser beams along the 
z-direction with a frequency difference around 100 kHz 
are applied to the atomic cloud and cause diffraction of 
the atoms. In the moving frame of the optical standing 
wave pattern, Bragg diffraction conserves kinetic energy 
and the atoms coherently populate two components at 
the incident momentum —hk (corresponding to zero mo- 
mentum in the laboratory frame) and at Hk, differing by 
twice the photon momentum. For suitably choosen pa- 
rameters the diffraction process is fast and interactions 
can be ignored. The diffraction is then a linear process 
and each atom is put in a superposition of remaining at 
rest, and having received twice the photon momentum. 
Allowed to propagate freely the atoms will coherently 
populate two spatially separated regions of space after 
less than a msec. If we neglect the role of interactions in 
this short phase, we can write the state of the system 




ni;(/)i(r))|n2;02(r)) (1) 



where |ni;0i(r)) denotes a state with rii atoms popu- 
lating the spatial wave funtion 4>i{r)\ n2 = N — m. In 
the subsequent dynamics, the wave functions (j)i{r) evolve 
with time, and we describe this evolution with the Gross- 
Pitaevskii equations 



(2) 



where hi is the single particle hamiltonian (the atoms 
may be free, in which case hi equals the kinetic energy 
operator, they may fall under gravity, or they may be 
trapped in a weak trapping potential) and g — ^nh^as/m 
is the collisional interaction strength (a^ is the s-wave 
scattering length). 

To compute the effect of collisions on the spatial dy- 
namics it is an adequate approximation to assume that 
half of the atoms N/2 populate each component. A more 



precise account of the interactions will be needed, how- 
ever, when we examine the details of the evolution across 
the binomial variation in n^. Effectively, the functions (pi 
and (j)2 define two modes for the atoms with creation and 
annihilation operators a\ and a^, and the dynamics asso- 
ciated with the distribution of atoms among the modes 
is accounted for by a two-mode Hamiltonian: 
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(3) 



where Ti = J (Pr\(j>i\'^, is determined from the solution of 
Eq.(||). It is not strictly correct to separate the spatial 
and the population dynamics, especially not if the con- 
densates are strongly interacting, which suggests a vari- 
ation of the preferred mode function (j)i over the binomi- 
nal distribution of rii around N/2. In the Thomas-Fermi 
regime, this variation has been taken into account within 
a broader ansatz for the state of a similar two-component 
system suggesting, however, that this effect is taken well 
care of if Ti is multiplied by a numerical factor of order 
unity M . 



B. Collective spin picture 

The Hamiltonian (|^) will cause a quadratic (with iii) 
evolution of the phase of the expansion coefficients in the 
state (0). This evolution is conveniently described by 
defining the collective spin operators 

S+ = 4ai, = a{a2, = (0^02 - a\ai)/2. (4) 

The Hamiltonian (^ can now be written: 

H = hxS', + f{N) (5) 

with hx = gT, where we have assumed for simplicity 
Ii = I2 = T. The f{N) term gives an overall phase to 
the system which we can neglect while the quadratic term 
5^ has been studied in detail by Kitagawa andUeda 0, 
who pointed out that precisely this Hamitonian leads to 
squeezing, i.e., reduction of the variance of a spin com- 
ponent orthogonal to the direction of the mean spin as 
illustrated in Fig. ^ 

Starting with the state (1), a spin coherent state 
with the macroscopic spin pointing along the x-axis, the 
Hamiltonian (5) produces a state which retains a large 
value of (Sx) and vanishing (Sy^z)- The initial variance 
Var(S'j,) = N/4 changes due to the interaction, and one 
identifies a particular component Se = cos 6 Sz — 
sinO Sy for which the variance is minimal. In terms of 
A = 1- cos^-H2xt) and B = 4sin(xi) cos^^^^^^) ^j^g 
angle 9 is given by 
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and the minimal variance is given by 
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FIG. 2: Phase- and population dynamics in the collective 
spin picture. The ^-component Sz represents the population 
difference between the two spatial components, and the ini- 
tial single component condensate starts out on the south pole 
with equal uncertainty in the Sx and 5*1, spin components. 
The preparation of the state with all atoms in a given su- 
perposition of the spatial modes 1 and 2 is represented by 
a simple rotation of the spin vector. The first Bragg pulse 
is a n/2 pulse and the spin is rotated to the equator of the 
Bloch sphere. By convention we choose the rotation axis to 
be the y-axis so the mean spin ends up along the a;-axis. The 
term in Eq. (^) now deforms the uncertainty elipse while 
leaving the mean spin along x. This deformation stretches the 
uncertainty elipse in the horisontal direction while Var(S'z) is 
kept constant. To obtain reduced fluctuations in the popula- 
tions of the two spatial modes, i.e. in Sz, a rotation by the 
angle shown in the figure is applied by a final Bragg pulse, 
chosen with approriate phases so that the rotation is around 
the direction of the mean spin. 



This spin squeezing is equivalent to a collapse of the 
relative phase, smeared out by the interaction |l^, 
It is important to notice that the phase collapse is not a 
decoherent process; the system is still in a pure state, and 
by a simple unitary operation described below (a spin ro- 
tation, cf. Fig. (H)) it can be transformed to a state with 
the desired reduced population fluctuations. We focus 
on population fluctuations, but it is worhwile mention- 
ing that it is equally possible to obtain a state where the 
relative phase of the two components has become well 
defined. 

The Hamiltonian (5) applies whenever the two compo- 
nents are spatially well separated, and we let the atomic 
clouds separate for a time interval of duration T (typi- 
cally some msecs). We then apply a new Bragg diffrac- 
tion pulse to induce a complete transfer between states 



with momenta ±?ifc in the moving frame, i.e. and 2hk 
in the laboratory frame, so that the two components now 
approach each other. This operation is carried out on the 
two spatially separated components, and it is accounted 
for by inclusion of the interaction with the Bragg fields in 
hi in Eq.(2). Note that with the convention that the col- 
lective spin refers to the spatially separated components, 
this Bragg pulse has no effect on the values of the spin 
components; the solutions (r, t) constitute a 'rotating 
frame' for our calculation of populations and coherences. 

After a second time interval of duration T with free 
evolution of the two condensate components, they again 
overlap in space. At this point it is possible to observe 
interference fringes in the spatial density profile plf . 
Note that although Fig. 2 is suggestive that spin measure- 
ments and rotations can be carried out at any time, it is 
only when the two spatial components overlap, the spin 
variables represent physically accessible quantities. At 
the time of full overlap we can apply a final beam split- 
ter Bragg pulse to recombine the two components into 
two new output components: 



%ragg = exp [ia (cos (j) Sx 



sin(f) Sy)] 



(8) 



Here we have written the Bragg evolution operator in the 
spin picture and as a 2 x 2 matrix giving the mode an- 
nihilation operators bi and 62 after the pulse in terms of 
the operators oi and 02 before the pulse. We are inter- 
ested in creating a state with two separate components 
with sub-binomial counting statistics, that is, we want to 
chose 4>, the phase of the pulse, and a, the duration, such 
that the counting statistics of 61 and 62 will benefit from 
the correlations in oi and 02. In the spin picture it is 
clear that this operation should be a rotation around the 
mean spin to align the squeezed direction with the z-axis 
which represents population differences. If we denote by 
Sa the spin component that the last Bragg pulse rotates 
into the z-direction then from the definition (||) we get 
the mean and variance of the populations of one of the 
output beams 



(ridot) - ^N+{Sa) 
Var(ndot) = Var(S'a). 



(9) 



If (/) = we rotate around the cc-axis and the angle a 
should ideally be choosen equal to 9 of Fig. ^ and Eq. (||). 
With no spin squeezing one finds Var(ridct) — N/A in 
agreement with the initial binomial distribution. 

From Eq. (||) it is clear that the distinction between 
different rotation axes, all lying in the xy-plane is equiv- 
alent to the one between complex phases (/> of the cou- 
pling amplitude, i.e. the relative phase of the two coun- 
terpropagating laser beams creating the moving stand- 
ing wave. Alternatively, we see that the location of the 
Bragg diffraction pattern, due to the translational prop- 
erties of momentum eigenstates, controls precisely the 
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phase of the couphng. In Ref. |T5|, we learned that the 
mean spin and the rotation angle are predicted with ad- 
equate precision in the simple model (||J|), so that appli- 
cation of Bragg pulses with these variables should suffice 
to yield substantial reduction of the population fluctua- 
tions of the two atomic outputs. It is also straightforward 
to optimize the parameters in experiments. The phases 
of the diffracting lasers are adjusted so that the two 
output beams have the same mean occupancy, indepen- 
dently of the duration of the coupling. The solid curve in 
Fig. ^ shows the squeezing factor (A^/4)/Var(ndct) (large 
for strong squeezing) as a function of a. With a spin 
squeezed sample the variance of the populations show 
a strong dependence on the duration of the last Bragg 
pulse, which should therefore be adjusted to identify the 
output with minimal fluctuations. 



III. EXPECTED RESULTS 

Let us turn to the discussion of the experimental fea- 
sibility of our proposal. Splitting and recombination of 
one-component condensates have been done in several 
laboratories, and the coherence properties have been ver- 
ified and explicitly utilized in a number of imaging exper- 
iments; our proposal follows closely the experiments at 
NIST, where the phase variation of the condensate gives 
rise to a density variation in the recombined output [ p2[ . 
These experiments have been done in the limit of inter- 
acting condensates which is of course essential for our 
proposal: In |^ the mean field repulsion was observed 
and it was shown not to prevent a nearly perfect overlap 
of the spatial modes at the recombination, and in |2|] 
a soliton was imprinted in one of the components to be 
subsequently detected in the output. We note, that also 
the phase variation around quantized vortices has been 
studied by similar interference imaging, both in the case 
where a vortex was prepared prior to the splitting and a 
dislocation appears in the interference fringes p3| , and in 
the case where a condensate was split and a vortex was 
subsequently created by stirring only one component [ p4[ . 

To assess the strength of the squeezing interaction we 
need to evaluate the parameter x of Eq.(||). In the 
Thomas-Fermi approximation, valid for an interaction 
dominated condensate, one finds the simple result: 



1= |</.rd^r = 



10 1 



TF 



(10) 



where Vtf = {4:T^/5)RxRyRz with Ri the Thomas-Fermi 
radius of the condensate in the i-direction. With the 
parameters of the experiments with condensates, being 
split and recombined by Bragg diffraction at NIST, [^ , 
one finds x — 4.7 x lO^^s^^ or a factor of ^ 4 re- 
duction in the population variances in the 2 msec du- 
ration of the experiment. This result is improved to 
a factor of more than 12 by increasing the duration 
to 4 msec. The theoretical maximal squeezing factor 



- (l/3)(iV/6)-2/3 ~ 1.3x10^ would require the com- 
ponents to be separated for ~ 6-^/^x^~^^^ ~ 200msec 
assuming a constant x throughout the experiment. In 
reality the condensate components expand if the trap is 
turned off during the experiment and the continuous re- 
duction of X over longer times according to Eq. (|l^) would 
have to be determined by an integration of the Gross- 
Pitaevskii equation (2). 



IV. EFFECT OF IMPERFECT OVERLAP 

As a final issue we wish to discuss the matching of the 
spatial wave functions of the two components. After the 
physical interactions described above, one will detect a 
number of atoms with momenta around zero, and a num- 
ber of atoms with momenta around 2hk, and the analysis 
shows that these numbers will fluctuate less than if they 
were given by a binomial distribution. This reduction is 
due to the correlations between the atoms following the 
lower and the upper paths in the interferometer, and it 
is hence important that atoms from the lower path, be- 
ing diffracted by the last Bragg pulse, occupy the same 
spatial state as the undiffracted component of the upper 
path, cf. Fig. |l|. In the extreme case where these two 
contributions are orthogonal and distinguishable, there 
will be no squeezing effect at all, and in the case where 
they overlap 100 %, the squeezing is given by the above 
simple analysis. 

Atoms leaving the interferometer at zero momentum 
can in general be described by a superposition of 02 j the 
undiffracted upper path wave function, and a state 02X 
orthogonal to (j)2 ■ This gives rise to two orthogonal modes 
which both contribute to the detected number of atoms 
around zero momentum. We can express the field an- 
nihilation and creation operators &2 , ^2J_ and b\ j ^Ll ^'^^ 
these two modes in terms of the operators ai,a2 and 
a\ , flj given above for the modes prior to the last Bragg 
pulse and operators for two additional 'vacuum' modes 
flu, ^21 and a\^, a^, which are necessary to ensure uni- 
tary of the beam splitter pulse. Letting c denote the 
overlap of the (normalized) diffracted lower path compo- 
nent with the upper path undiffracted component and s 
the overlap with the orthogonal complement 



— 2ikx 



01 (r) = c(j)2(r) + s(j)2L{r) 



(11) 



(|sp +J,cP — 1)- In terms of a and of the last Bragg 
pulse (H) we have: 

62 — cos — a2 + i sin ^e"*"^ (c ai -f s an) (12) 
621 = cos — 021 + i sin — e^*"^ (s* ai — c* a-]±) . (13) 

Focusing on 62 we see that for |c| < 1 we can still control 
the mixing and pick the squeezed combination of Oi and 
a2. However, some contributions from the unoccupied 
modes have been introduced, and the situation is similar 
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to the detection of squeezed light: A non-unit overlap v 
between the squeezed mode of light and a detector mode 
as, e.g., modelled by a beam splitter transmitting a frac- 
tion ji^p of the squeezed light to the detector, leads to 
admixture with a vacuum field and the deterioration of 
the squeezing, Var(X)dot = \v\^N^Y{X)sq + (1 - ji'P), 
where the X is the scaled quadrature field component 
with unit variance in the vacuum state. 

In a real experiment it is difficult to count the atoms in 
mode (^2 exclusively. It is more realistic that one will only 
have access to the total number of atoms in 02 and (/)2j_ . 
We then have to include 621 in the analysis, and if the 
phase of the Bragg diffraction is chosen to compensate 
the phase of the complex overlap c, we find: 



{blb2 + bl^b2L) 



(cos a Sz 



IcI sin a Sy 



N 



N 
1' 



(14) 



in terms of the spin component Sa' 'prior to the last 
Bragg pulse. The angle o! is related to a, the rotation 
angle given by the duration of the last Bragg pulse, by 
tana' = |c| tana, and wc have introduced the factor A = 

\l |cp/(sin^ a' -f |cP cos^ a'). 

The spin component Sa' = cos a'5'; 
direction perpendicular to the mean spin. We therefore 
have {Sa') = and (ndct) = N/2 independently of a 
The variance of Tidot is given by 



sin a' Sy is in a 



Var(ndct) = ("dot) - ('Met) 



N 



= A"Var(5„0 + (l- A')-. (15) 

If the overlap is complete |c| = 1 = A, we can pick out the 
maximally squeezed spin component to yield the minimal 
variance of ndct- In case of a non-perfect overlap, we 
see that Eq. (15) introduces a non-squeezed contribution, 
scaling as N/T, and it makes the number fluctuations 
depend on the noise of the spin component Sa' rather 
than of Sa ■ We therefore want a' to be close to the actual 
direction of squeezing 9 in Fig. 2 and Eq.(||), which is 
obtained by choosing a Bragg pulse with longer duration 
than what is optimal in the ideal case. 

In Fig. ^ we show the variation of (7V/4)/Var(7idct) 
with a for different values of the overlap c. When |c| 
is reduced the squeezing factor drops, but the curves il- 
lustrate that the results improves if one choses a larger 
rotation angle a. 

Note that the mode function mismatch does not have 
the same detrimental effect as in the detection of a 
single mode squeezed field: the admixture of the vac- 
uum modes, and hence the vacuum contribution N/A 



to Var(ndot), does not scale with \s\ 



1 



but 



rather with \s\^ sin^ a, and for small rotation angles this 
is a small number. For strong squeezing the uncertainty 
elipse is almost horisontal and only rotation by a small 




FIG. 3: The reduction of Var(7idGt) is illustrated by plotting 
the reciprocal ratio (A''/4)/Var(ndot) as a function of the aj 
plied Bragg pulse rotation angle for parameters as in Sec 
with 4 msec separation. The results are shown for different 
values of the overlap, c = 1 (solid line), c — 0.9 (dashed line), 
c = 0.8 (dash-dotted curve), and c — 0.6 (dotted curve). The 
value of a at the maximum of the c = 1 curve displays the ori- 
entation of the squeezing elipse, i.e. 6 of Fig. |^. As is apparent 
from the other curves, one may compensate for a non-perfect 
overlap by applying a larger Bragg pulse rotation. 



a is needed to observe squeezing in the number fluctu- 
ations. In this case even a significantly reduced overlap 
does not prevent a significant noise reduction for ridot- 
This is illustrated when we write the maximal squeezing 
factor attainable for given c and given Y0.r{Se)- 



Var(ndct) = Var(5e) 1 + (1 - jcp)- 



f - Var(^g) 
c|2f +Var(5e) 




. ^ (Var(5.)«^). (16) 

Iq) we assume that the spin state described by 
is a minimum uncertainty state). 

V. CONCLUSION 

We have proposed to use the coUisional interactions 
between atoms in the separate arms of a spatial inter- 
ferometer to suitably manipulate the amplitudes on the 
different number state components. We analyzed the 
achievements of current condensate interferometry exper- 
iments, and we found that significant number squeezing 
should be found in the interferometer output, i.e., the in- 
terferometer acts as a quantum beam splitter for atoms 
with unique noise properties. Good overlap of the out- 
put mode functions is necessary to observe interference 
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effects in the interferometer and it is also important for sizable noise reduction is possible with interferometeric 
the noise reduction in our beam splitter. We showed, overlaps easily achievable in experiments, 
however, that the overlap is not a critical parameter and 
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